Electrostrictive fluid pressure from a laser beam 
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Recent times have seen surge of research activity on systems combining fluid mechanics and electromagnetic 
fields. In radiation optics, whenever information about the distribution of pressure in a dielectric fluid is 
required, the contribution from electrostriction becomes important. In the present paper we calculate how 
the local pressure varies with position and time when a laser beam is imposed in a uniform fluid. A Gaussian 
intensity profile of arbitrary time dependence is assumed for the beam, and general results are derived in 
this case. For demonstration we analyze two different cases: first, that the beam is imposed suddenly 
(mathematically in the form of a step function); secondly, that the beam is switched on in a soft way. In both 
cases, simple analytical expressions for the pressure distribution are found. 
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I. INTRODUCTION 

Over the last decade, considerable experimental re- 
search has been conducted into the interplay of optics and 
fluid mechanics. The field of optofluidics, the application 
of microfluidic flows for optical purposes, is a topic in its 
infancy and is seen by many as holding very considerable 
technological promise^. Appreciable progress in optical 
manipulation of microflows has also been achieved; some 
examples are found in Refs. EHZl These impressive exper- 
imental developments create the need also for improved 
theoretical understanding of liquid-laser interactions. As 
a step towards this end we have investigated herein the 
time-dependent electrostrictive pressure due to a laser 
beam propagating through a liquid. 

The electrostriction force density is a local force den- 
sity acting within all dielectric fluids and solids subjected 
to an electric field, proportional to gradient of the field 
intensity. The force acts towards higher field strength, 
thus causing a local pressure increase in the center of a 
laser beam propagating through the medium. From a 
general viewpoint, when including also solid elastic me- 
dia, the electrostrictive susceptibility is strongly related 
to the elasto-optic susceptibility as they both represent 
the interaction of electric fields and an elastic deforma- 
tion in the medium. A very general account of media of 
this kind can be found in Nelson's boolP, and a compre- 
hensive statistical mechanical theory is provided by Rasa- 
iah, Stell and co-workerJ^R When restricting ourselves 
to electrostrictive phenomena in fluids we may mention, 
for instance, the early paper of Brueckner and JorngP^l 
in which it is discussed how the electrostrictive coupling 
between laser light and the medium can be important 
for rapidly growing instabilities, and hence may lead to 
large density fluctuations. Another work dealing with 
electrostriction in fluids is the series of papers by Zan- 
chini and Barett a 13 " 15 ^ here a discussion is given of the 
role of electrostriction for the attractive force between ca- 
pacitor plates completely immersed in a dielectric fluid. 

A work of particular interest in our context is the mi- 
crogravity experiment of Zimmerli et at on electrostric- 
tion in a fluid under extraterrestrial conditions^. This 



work is actually a generalization of the classic elec- 
trostriction experiment of Hakim and HighairP^ (The 
Hakim-Higham experiment is discussed more closely in 
Ref. [T51 section 2.3.). The last-mentioned authors mea- 
sured the pressure between closely spaced electrodes im- 
mersed in a nonpolar liquid (carbon tetrachloride and n- 
hexane). The essential new element in the experiment of 
Zimmerli et al. was to examine electrostriction in a liquid, 
sulfur hexafluoride (SFg) near the liquid-vapor critical 
point, T c = 319 K, where the compressibility was seven 
orders of magnitude higher than in the liquids studied by 
Hakim and Higham. Because of this compressibility, even 
a moderate electric field produced an easily measurable 
density change in the nonpolar fluid. To avoid distur- 
bances from stratification gravity effects on Earth, the 
experiment was carried out in outer space, in the Space 
Shuttle Columbia in 1994. For the application of electro- 
dynamic theory to fluids, the experiment is important. 

However, generally speaking the occurrence of an elec- 
trostrictive pressure term in a dielectric fluid has for the 
most part received only moderate attention in practi- 
cal situations. The reason for this is quite evident: the 
electrostrictive force density is generally written as a gra- 
dient, and the total electrostrictive force on a dielectric 
specimen, which may be written as an exterior surface 
integral, gives no contribution to the force on the spec- 
imen as a whole if the surroundings are non-dielectric. 
Under usual physical conditions, the electrostrictive ef- 
fect therefore only becomes of importance in considera- 
tions of local pressure distributions in matter, accessible 
experimentally usually via optical methods. 

Although deformations and flows generated by elec- 
tromagnetic fields can generally be described without in- 
clusion of electrostriction, this may not be so for their 
stability. This point was emphasized by one of us al- 
ready some time ageff^ and we return to it in the follow- 
ing few examples. The deformations of fluids by lasers, 
under static conditions as well as under time-dependent 
conditions, are therefore adequately described in terms of 
Maxwell's stress tensor, without electrostriction included 
at all. See, for instance, some of the Bordeaux papers in 
Refs. |H and [7J and also some related theoretical pa- 
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pers in Refs. |2"U1 and |2"T1 For full stability considerations, 
however, electrostriction could play an important role. 

Electrostriction can moreover be made useful, for ex- 
ample in laser-induced thermal acoustics (LITA) which 
has been made a tool for measuring a ra nge of quantities 
including the speed of sound in a mediuirPSGS] The elec- 
trostrictive force, being able to effectuate local changes in 
density, can also locally alter a medium's refractive index, 
giving rise to non-linear optical response even at moder- 
ate field strengtfP21, resulting in such effects as anomalous 
Stokes gairP2l and optical self-focusing^. 

A few examples may serve to illustrate the typical role 
of electrostriction. Firstly, it is noteworthy that the elec- 
trostrictive force is of the same order of magnitude as the 
usual electromagnetic force density — (e /2)£ ,2 V£ acting 
in regions where the permittivity e varies with position, 
typically in dielectric boundary layers. (We write the 
constitutive relations as D = £o£E, B = HofiH., so that e 
and n become nondimcnsional.) Consider for definitcness 
the usual textbook situation in which two metal plates 
are partially immersed in a dielectric liquid (the situa- 
tion is discussed, e.g. in section 2.2 of Ref. IT8"|) . When 
there is an electric field between the plates, the liquid is 
known to rise to a height h that can be easily calculated. 
The main physical principle here is the balance between 
the lifting surface force at the free surface and the grav- 
itational force. A nontrivial point is, however, the effect 
of the electrostrictive force: its presence is necessary to 
secure that the excess electric pressure in the liquid is 
strong enough to make it possible for the liquid to rise 
as a coherent whole 19 . The electrostrictive effect thus 
serves to stabilize the system - it is intimately related to 
thermodynamic stability - but it does not play a role for 
the magnitude of the height h. 

The observability of electrostrictive phenomena de- 
pends on one single physical variable in the system, 
namely the velocity of sound. In the classic experiment 
of Ashkin and Dziedzic^, our second example, a narrow 
beam of light was sent through a dielectric liquid and an 
elevation of the free surface (in the order of 1 p.m) was 
observed. The mean time between molecular collisions in 
such a liquid is of order 1 ps, and sound needs only a time 
of about 7 ns to traverse the cross section of the beam 
(waist radius wq = 4.5 /mi). This means that the elastic 
pressure has had sufficient time to build itself up long be- 
fore experimental effects become visible. As shown from 
a detailed calculation in Fig. 9 in Ref. HH1 the maximum 
dynamic electrostrictive pressure occurred at times t ^3 
ns after the sudden onset of the pulse at t = and died 
out after 8-10 ns. By contrast, the hydrodynamical mo- 
tion of the free surface needed hundreds of nanoseconds 
to develop. For practical purposes it is thus legitimate 
to introduce the concept of a hydrodynamical pressure 
p, and to assume that the electromagnetic forces start to 
act simultaneously throughout the liquid. 

A third example of essentially the same kind is fur- 
nished by the radiation pressure experiment of Zhang and 
Chang from 1988^1 the authors observed the radiation- 



induced deformation of the surface of a micrometer- 
droplet illuminated by a laser beam. Again, the elec- 
trostrictive pressure in the interior of the droplet is coun- 
terbalanced by an elastic pressure having had sufficient 
time to build itself up. The deformation of the droplet 
can thus be calculated as though there were no elec- 
trostrictive effects at all. The theory of this experiment 
has been worked out in Refs. and HH) 

Under special circumstances it becomes possible to 
make the electrostrictive effect discernible in force mea- 
surements even without relying upon optical meth- 
ods. The Goetz-Zahn experiment belongs to this 
category^2H22| ^he main principle of this experiment was 
to apply a high frequency electric field E = E cos uit be- 
tween the two plates of a condenser filled with a liquid 
(polar or nonpolar) , and to measure the appropriate har- 
monic component of the force between the plates. Strong 
fields were necessary for this purpose, of the order of 
E = I0 6 V/m. The reason why the electrostrictive part 
of the force was detected, was that the liquid was in a 
state of mechanical non-equilibrium; the fields were os- 
cillating so quickly that the elastic counter pressure had 
insufficient time to build. To propagate a distance of 1 
cm, sound needs a time of about 10 |_is, which is more 
than the period 2ir/us ~ I /xs of the oscillations. This 
case was discussed in detailed in Ref. [T51 section 2.4. 

In the following we will calculate the distribution of 
the electrostrictive pressure in a liquid under typical ex- 
perimental conditions. The practical utility of such an 
undertaking lies in the increasing need to have detailed 
information about the stress distribution in soft matter 
in the presence of lasers. In view of the increasing in- 
terest in liquid-laser interactions at a microscopically de- 
tailed level, such a calculation is timely. We consider for 
demonstration two different time evolutions of the laser 
beam; the case in which the pulse is switched on and off 
abruptly, and the case of a softly peaked pulse with expo- 
nential tail, while our final result is valid for an arbitrary 
time evolution and can be of use, for example, in ap- 
plications within laser-induced acoustics. In both cases 
relatively simple expressions are obtained for the pres- 
sure distribution in the liquid as a function of distance 
from the beam axis and time [corrections due to beam 
divergence are ignored] . A natural time and length scale 
for the problem are provided by the beam width and the 
velocity of sound. Our final results in the two cases are 



Eqs. (33) and (51), in which the pressure distribution is 



given as function of time and radius rescaled with respect 
to these scales. 



A. Basic formalism 

We regard a laser beam propagating along the z di- 
rection in a homogeneous liquid, with beam intensity de- 
scribed by a Gaussian so that the squared electric field 
in cylindrical coordinates r = (p, 6, z), 



E 2 (r,t)=E 2 (v)T(t) 



(1) 
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FIG. 1. The situation considered: A Gaussian laser beam 
of width wo > A propagates through a medium of relative 
permittivity e (A: wavelength). The laser intensity varies in 
time according to a function T(t). 



with (c.f. e.g. Ref. [18] and chapter 16 of Ref. 



E 2 (r) = 



2P 



exp 



2p 2 

w 2 {z) 



(2) 



irsoncw 2 (z) 

The situation is shown in figure [T] Here, P is the total 
incident laser power, n is the liquid's index of refraction, 
and because of inevitable beam divergence its radius is z 
dependent and of the form 



w(z) = wq\J 1 + z 2 /l 



(3) 



where Wq is the waist radius and In = irw 2 /\ is the 
Rayleigh length with A the wavelength of the laser light 
in the liquid, and c is the velocity of light in vacuum. 
In the following we shall assume the beam width wq to 
be much greater than the laser wavelength A, so that 
w(z) sa wq and remains approximately constant over a 
propagation distance of many beam diameters. 

The electrostrictive force due to the presence of the 
beam is V% where 



X(p,t) = \e E 2 {p,t)g r 



An 2 
dg„ 



X (0)T(t)exp 



2p^ 
w 2 (z) 



where g m is the fluid's density and x(0) is the initial and 
original value of X- The speed of light is much greater 
than that of sound, hence we can assume all changes in 
the electromagnetic laser field to be instantaneous and 
simultaneous throughout the liquid. We then clearly have 
X(p, t) — x(p)T(t). Assuming the liquid to obey the one- 
component Lorentz-Lorenz (Clausius-Mossotti) relation 
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e n 2 + 2 
where a is molecular polarizability gives 



X (p) = \e E 2 {p){n 2 -l){n 2 + 2). 



Explicitly, 



x(o) = 



(n 2 - l)(n 2 + 2)P 



3irncwQ 



(5) 



(6) 



(7) 



It is to be noted that the Clausius-Mossotti relation, ac- 
cording to which the permittivity depends on the fluid 
density only, refers to a non-polar fluid. The relation is 
convenient to use in practice, as one does not have to dis- 
tinguish between the isothermal derivative (dn 2 /dg m )T 
and the adiabatic derivative (dn 2 /dg m )s- It is found to 
apply with good accuracy for gases, for many liquids, 
and even for some solids. For polar liquids, in order to 
get better accuracy, the Onsager relation is available^ 
with various modifications^. In our theory the use of 
another relation for n(g m ) amounts only to a different 
form of x(0), and since we will work to linear order in 
the density variations in the liquid, x(0) will enter only 
as a prefactor throughout the below theory. Hence the 
choice of relation n(g m ) does not affect the following cal- 
culations and results in any essential way. 

We ought to point out at this point that we ignore 
other effects besides electrostriction, in particular, the 
quadratic electro-optic effect. It is known that in strong 
electric fields a liquid can become optically anisotropic, 
the polarizability ay parallel to the field being different 
from the polarizability cvj_ in the transverse direction. 
The field dependence of the polarizability in turn influ- 
ences the refractive index through the Lorentz-Lorenz re- 
lation (the optical analogue of the Clausius-Mossotti re- 
lation). On may estimate the magnitude of the quadratic 
electro-optic effect by considering the nonlinear quantum 
theory of the refractive index in an isotropic medium. 
For instance, regarding an artificial medium consisting 
of two-level atoms (see p. 147 of Ref. [TBI and further ref- 
erences therein), it follows that an = «o/(l + E 2 /E 2 ). 
Here ao is the polarizability in the absence of the field 
and Eq is the constant Eq = huj s j (y/2exo s ), where Tiuj s is 
the energy difference between the two energy levels of the 
atom and xq s the matrix element of the position opera- 
tor between the same levels. For simple non-polar liquids, 
typically w s ~ 10 s , so that if the matrix element is 
in the order of the Bohr radius, we obtain E ~ 10 10 to 
10 11 V/m. This is a very high value, implying that the 
quadratic electro-optic effect can be safely ignored under 
usual physical conditions. 

To study the pressure that builds due to the elec- 
trostrictive force, consider the linearized Euler equation 



g m d t v = -Vp + Vx 



(8) 



[dt = d/dt]. It follows from Q and the fact that the 
same equation must hold before the arrival of the beam, 
that V x v = 0, and hence we may write the velocity field 
in the form of a velocity potential, 



(9) 



Thus, because of gauge invariance of potentials we get 
the simple equation 



p(p, t) = -g m d t <Z>(p, t) + x(p, t). 



(10) 



Writing p(p,t) = p a + p'(p,t) and g m (p,t) = g m0 + 
g' m (r,t) where po,g m o are the constant, initial pressure 
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and density, and p' <C Po, g' m ^ Qmo, we linearize w.r.t. 
p', g' m and go to the gauge frame, 



P'(P,t) = -Qniod t $(p,t) +x(p,t). 



(11) 



The second term on the right hand side is due to the 
direct electrostrictive force from the local electric field 
density whereas the first term comes from propagation 
of pressure waves from the surrounding area. Spherical 
pressure waves are emitted from all points within the 
beam when the electric field varies with time. [Field 
variations over an optical period is of no consequence; 
the material response time is much longer than this, and 
these periodic field variations simply average to zero.] 
With the continuity equation 



d t Q m + V • (g m v) = 



and noting that 



P' = u 2 g' m 



(12) 



(13) 



where u = [{dp' /dg^s] 1 ^ 2 is the speed of sound, the 
governing equation for <&{p,t) is obtained 



(V 2 - u- 2 d 2 Mp, t) = fdtxip, t). (14) 



II. STEP-FUNCTION SWITCH-ON OF BEAM 

Assume first that the beam is switched on suddenly at 
t = and off again at to. 



T(t) = Q(t)-e(t-t ) 



(15) 



where is the unit step function. Because the velocity 
of light is immensely greater than the velocity of sound, 
we may safely assume the onset of electrostrictive forces 
to be simultaneous throughout the liquid. We have 

dtx(p,t) = x(p)d t T(t) = X (p)[8(t) - 5(t - to)]. (16) 



We define the Green's function of Eq. (|14p , which by 
definition satisfies 

(V 2 - u- 2 d 2 )G{v, r', t, t') = S 3 {r' - r)6(t' - t) (17) 

which has solution (e.g. Ref. [31)1 p. 244) 

1 



G(v',r,t,t') 



4?r r' - r| 



6 [t 1 - 1 + 



(18) 



This solution requires that scattering of sound waves off 
fluid surfaces may be neglected. In microscopic systems 
where such effects are of importance, a Green's func- 
tion taking these boundaries into account is required - 
Green's functions for the scalar wave equation have been 
found for a wealth of geometries. 



From (14 1 and (18) it follows that (noting 5(ax) 
6(x)/\a\) 
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g m ou 
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At' / d\'G{v',v,t,t')d tX {r',t') (19) 



(/dpi 



/e xp(-2p' 2 / W 2 ) 



R 



x [S(tu - R) - 5{ut - ut - R)] 
where we define 

R=|r'-r|; R = |R| 



(20) 



(21) 



Since the geometry is obviously symmetric under z — > — z 
and 9 — > —9 we have halved the integration ranges of z' 
and 9' above. For the same reason we are free to choose 
r to lie on the x axis without loss of generality: z = 0, 
9 = 0. To integrate out the remaining delta function we 
integrate with respect to R using 



R 2 = P 2 



P' 



to substitute 



dz' 
~R~ 



dR 



2pp' cos 9' 



dR 



y/R 2 - p 2 - p> 2 + 2pp>cos9> 



(22) 



(23) 



The delta functions in Eq. (20) pick out two spherical 



shells of values of (r',9') centered at r and having radii 
ut and, for t > to, u(t — to), respectively. Thus, consider- 
ing the term of radius ut [the second term is essentially 
identical, as is clear from Eq. (20)] 



d9' 

n JO 

d9' 



d9' 



p'dp'f(p') 



p'dp'f(p') 



dz' 
~R 



6(tu - R) 



dR 5(R - tu) 



P'dp'f(p') 



Vu 2 t 2 - TV 



e(ut-K) 



where we define 

n 2 {t) = P 2 



- 2pp'cos9'. 



(24) 



(25) 



In the above, the exponential in ( 20 1 depending only on 
p' were gathered in a shorthand function f(p'). 

The function restricts the integration area to a cir- 
cular disc of radius ut centered at (p, 0, 0). It is natural to 
integrate instead by parametrising according to a cylidri- 
cal coordinate system where the center of this disc lies at 
the origin as shown in figure [2] We note that m a new 
coordinate system (1Z, ip), 



TZ 2 + 2pK cos (p. 



(26) 



Now the G function is equivalent with the upper integra- 
tion limits, and can be removed. An equivalent expres- 
sion is obtained for the second term, but with t — > t — to, 
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expressions 
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FIG. 2. Coordinate shift: translating the z axis to the point 

(p,o). 



and only allowing contributions for t > t for reasons of 
causality. 

Using the relation 



&P e -aco S <p = / ( a ) 
7T 



(27) 



where ioM is the modified Bessel function of the first 
kind of order zero, we thus obtain a general expression 
for $(p, t): 

wm f /of^lexpf- 2 ^ 7 ^ 
Mo t ) = A J em / t?h7? iv V 



e(t-to) / ^d^ v v 



(28) 



The gauge pressure can now be calculated Eq. ( 11 ), and 



we shall do so in a moment. First, we will simplify our 
notation considerably by introducing rescaled variables 
for time and length. 



A. Rescaled time and radius 

For further analysis we introduce the typical time and 
length scales 



- - ^° ■ i - W ° 



(29) 



satisfying ut = p. The corresponding dimensionless time 
and space variables we define as 



e = p'/p; t = p/r, 

t = t/t; t = to/i; 



(30a) 
(30b) 



and substituting for 1Z — > x = x = \Ju 2 t 2 — 7Z 2 /p as 
well as rescaling we get the much simpler looking general 



z X(0)i 

Qrntt 



$(£,t)-$(£,t-t ) 



$K,t) =6(r)e^ 2 f dxI (2^Vr 2 
Jo 



x 2 )e x 



(31a) 
(31b) 



The integral in Eq (31b) is simple to evaluate numeri- 



cally. With the re-scaled variables and the dimension- 
less $ function defined above, Eq. (11) becomes [the 
re-scaled time envelope function is obviously T(t) = 
6(r) - 0(r - r )] 

Pit r)=- x(0)[d T $(£, t) - a r $(e, t - to) 

-T(T)e-« 2 ]. (32) 

Evaluating the derivatives we finally find the pressure 
distribution as function of (rescaled) radius and time 



P(t>r) =x(0)[P(£,t)-P(£,t-t )] 
rT xdxe~ 



P(£,r) =2Te" 4 0(t) 



o \/t 2 — X 2 



(33a) 



(33b) 



where in the last line we have substituted x — > y/r 2 — x 2 . 



Equations ( 33 1 are the final result of the pressure calcu- 



lation in the case of step-function switch-on and -off of 
the laser. 

The generalisation to a series of N pulses turned on at 
times ti and off at times £rjj (U < t(si < U+i) is obvious: 



P(£, r) = X (0) J2l P & ^ - T - 



(34) 



A plot of the functions r) and P(£,t) are shown 
in figure [3j One notes how the pressure on the axis rises 
steeply when the laser is switched on and then relaxes 
monotonously towards a new equilibrium value which has 
the very simple form 



P{i,T) > C 5 



(35) 



We show this in the appendix. While it is not entirely 
trivial to obtain the limit from the more general ex- 
pression, this result could have been obtained in a very 
straightforward way by assuming a static laser beam 
from the beginning. After reaching the peak value, the 
pressure relaxes quickly and the stationary situation is 
reached already at around 3t. The left panel of figure [5] 
also shows how a pressure wave travels outwards from the 
cylinder axis after switch-on at velocity u, and similarly, 
a wave of negative gauge pressure at switch-off. This fact 
allows closely controlled acoustic waves to be created by 
a laser beam, the principle behind laser-induced t herm al 
acoustics, already used for measurement purpose d 22 ! 23 ^. 
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FIG. 3. Left: the function r) for to = 5. Right: The reduced pressure p(£,i~)/x(0) for tq — 5 as given in Eq. (33b I 



B. Pressure time variations on the axis 



III. LASER PULSE WITH GENERAL TIME EVOLUTION 



A simple analytical expression can be found for points 
on the z axis, £ = 0. Then we are left with 



$(0, r) = 9(t) / dx exp(x 2 - r 2 ) = 9(t)F(t) (36) 
where F(x) = e~ x2 dse s2 is Dawson's integral. Insert- 



ing this into Eq. Q with F'(x) = 1 - 2a;F(a;) we get 
simply, in terms of the dimcnsionless units Eq. ( 30 ) , 



$(0, t) =Q{t)F(t) - 6(r - t )F(t - r ); (37) 
p(0,r) =2 X (0)[6(t)tF(t) 

-e(T-T )(T-T )F(r-T )]. (38) 

The dimensionless pressure p(0, r) /x(0) is shown in figure 
|4j The special case of the cylinder axis was considered 
also previouslj^, in agreement with Eq. (pressureAxis) . 




Whereas we considered a sharp turn-on and turn-off 
of the beam in the previous section, here we will be con- 
cerned with a more typical pulse transient. Within the 
same formalism as employed above, we are able to keep 
the function T(t) general throughout, yielding results of 
more general validity. For numerical examples in the fol- 
lowing, we choose the following typical form for a soft 
laser pulse, 



e 2 r 2 



(39) 



so that to = tTQ is again the duration of the pulse. Here 
and forthwith we make use of the rescaled variables de- 



fined in Eq. 30 We have 



T'(r) 



4r V r o 



-t/to, 



T"(r)=^l2-4 



-t/to 



(40) 
(41) 



The velocity potential is again given by Eq. (19), where- 
with we get after rescaling 



1 

47T 



-€' 2 



dr' / d£' S(t' -t + s)T'(t') 



1 1 d£'- T'(r-s) 



47T 



(42) 



FIG. 4. The pressure p(0, r) on the axis, Eq. q38J) , as function 
of dimensionless time r, divided by the prefactor x(0)- The 
pulse is turned on at r = and off at r = ro = 5. 



where we define 



s = £' - £; s = |sj 



(43) 



As before, we are free to choose £ = (£, 0, 0), and note 
that the integrand is invariant under rotation about the 
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£4 axis. The integral over £' runs over all of space, so 
we choose integration coordinates instead to be spherical 
coordinates with the origin at £ in the system and with 
the £' x axis as polar axis with polar angle 6 S . In the new 
coordinate system (s,8 s , ip s ) the integral reads 

*te> T ) = 4^J sdsT'{T-s)J^ d Vs J sm9 s dd s 
x exp[-£ 2 - s 2 (l - sin 2 S sin 2 - 2£scos6> s ] (44) 

where the exponent is — £' 2 (£, s, 6 S , (f> B ), as can be seen 
geometrically from figure [5] when noting that the length 
z' can be expressed as z' = ssmQ s sm<f> s and that |£'| 2 = 
£ 2 + .s 2 + 2£scos6» s = £' 2 + z' 2 . 



for the step-function pulse thus serve as derivation of the 
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FIG. 5. Coordinate transformation to spherical coordinates 
with origin at (£, 0, 0) in the cylindrical (£', 6 1 ', 2') system, and 
polar axis x'. 



We now make use of the relation 

d</) s e s2 si " 2 e = sin2 = 27re^ 2 sin2 e =/ (i S 2 sin 2 6.) 

Jo 
giving 



dsT'(r- s)e~^ 



x / d^e-^/olK* 2 -/9 2 )] cosh 2^/3 (45) 



*(£,t) 



where we substituted (3 = s cos # s . 

As a consistency check we should be able to insert 



the step-function time evolution derivative of Eq. (15 1, 
T'{t ) = S(t) — 8(t — t ) and get back the expression 
Eq. (31b). One quickly sees that this is so, provided 



dxe x ' I (2^y/r 2 -x 2 ) 



dpe^-^IvWiT 2 -/3 2 )] cosh 2^/3 (46) 



for arbitrary t and £. We have verified the identity nu- 
merically in the positive-positive (r, £) plane. The two 
different derivations of the re-scaled velocity potential 



potentially useful relation Eq. ( 46 ) . 

As in the previous case, the expression can be simpli- 
fied on the axis. Let us first use Eq. (46 ) to write Eq. (45 ) 

as 



-e 



dsT (t — s)e 



f dxe x2 h{2^s 2 ~x 2 ) 1 (47) 
J 



an expression which is still valid for an arbitrary T(r). 
On the z axis (£ = 0) $ becomes simply 

$(0,r)= / dsT'(r- s)F(s) (48) 
Jo 

where F(x) is again Dawson's integral. The functions 
$(0,t) and P(0, r) are plotted in figure [Hj 



<5(0,r) 

P(0,t) 




FIG. 6. The functions $ and P on the symmetry axis for 

T() = 5. 



The expression for the pressure is simple to write down 



from the generalization of Eq. (32 1, 



pUr)=x(0)[T(T)e-« -0 r *(£, T )]. (49) 



Analogous to and for comparison with Eq. ( 33a I , we de- 
fine 



K£,t) = x(0)P(£,t) 



(50) 



now with 

P(£,rHe-« 2 



T(r) 



f dsT"{T- s)e~ s2 f dxe x2 I (2^ 
Jo Jo 

T'(0)e- r2 f dxe x2 I (2^Vt 2 ~ x 2 ) 
Jo 



X 2 ) 



(51) 



This the expression for the pressure given a general T(r), 
although this form clearly requires the existence of T"(0) 
which is not satisfied for the step function considered 
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in section [nj For the soft pulse transient of Eq. ( 39 1 , 
T'(0) = and the last term vanishes. 

We plot $(£, r) and P(£, r) using the pulse transient 



medium, a principle already used in measurement set- 



( 39 1 in figure [7] Notice how much longer it takes for 
the gauge pressure to relax to zero again compared to 
the step function turn-off shown in figure [3] The be- 
haviour is otherwise comparable to that of figure [3] The 
generality of Eq. (51) implies for example that one can 



derive the exact acoustic pressure wave field throughout 
the medium which is set up by a train of laser pulses with 
known time dependence [reflection/absorption at mate- 
rial boundaries not included, of course - this would re- 
quire a different, geometry specific Green's function in 
Eq@]. 

Also in this case the expression takes a simpler form 
on the symmetry axis, where we are left with 



P(0, r) = T(t) - T'(0)F(t) - f dsT"(T - s)F(s) 

Jo 

where F(s) is again the Dawson integral. 



(52) 



IV. CONCLUSIONS 

We have calculated the local gauge pressure distribu- 
tion due to electrostrictive forces as a function of time 
and position, in the presence of a Gaussian laser beam 
propagating through a bulk fluid. The resulting expres- 
sions are simple and may form a natural part of an ana- 
lytical toolkit for experimental interpretation. A typical 
scale for length (p) and time (t) emerge, given by the 
beam width and the velocity of sound. Our final results 
are valid for a general time evolusion of the propagating 
laser beam, and for demonstratino purposes we consider 
two typical cases: a sudden switch-on/switch-off of the 
beam power, as well as a soft laser pulse. Reasonably 
simple expressions are obtained for both cases, as given, 
respectively, by our Eqs. (33 1 and (51) in terms of the 
re-scaled distance from the beam axis (£ = p/p), and 
time (t = t/t). On the beam axis, these expressions may 
be simplified further. The pressure relaxation time, the 
duration before the gauge pressure has returned to zero 
as the laser pulse is turned off, is surprisingly different in 
the two cases for the same pulse duration. 

Provided a stationary laser is switched on from zero 
to full intensity in a time shorter than the typical time 
scale t = wo/\/2u (wo'- beam waist, u: speed of sound), 
our result for a step function switch on are valid. After a 
transient period in which the pressure rises sharply and 
then relaxes to a stationary value, equilibrium is reached 
after a time of 3 — At, roughly the time for a sound wave 
to propagate the beam diameter. Under most circum- 
stances fluid motion is much slower than this, and elec- 
trostrictive pressure can be considered to be set up in- 
stantaneously. A cylindrical acoustic pressure wave trav- 
els outwards from the beam on switch-on, and similarly 
a wave of lower pressure is emitted at switch-off, allow- 
ing for highly controlled laser-induced acoustics in the 



22 23 



up d^ ' ^"- 1 '. Our result (51) provides in a simple way the 
acoustic pressure field which results from an arbitrary 
train of laser pulses anywhere in the medium at any time. 

Optofluidics and electrohydrodynamics are fields of re- 
search seeing impressive advances and promises many ap- 
plications in optics, sensing and measurement, microflu- 
idics and microchemistry, see e.g. Refs. [T] and [5] and ref- 
erences therein. So far this field of research has been 
driven by impressive experimental advances with theo- 
rists working to explain the observed phenomena. In this 
paper we provide another theoretical building block to- 
wards this end. Electrostriction, while not typically re- 
quired to describe electromagnetically induced fluid de- 
formations and flows themselves, is often of importance 
to their stability. This is examplified by perhaps the sim- 
plest example of electrohydrodynamics, that of a liquid 
rising between condensator plates due to the electric field. 
While the elevation can be calculated without reference 
to electrostriction, the latter plays the tacit but crucial 
role of ena bling the liquid column to rise as a coherent 
whole^ED. ^he methods and results presented herein 
could have many applications in analyses of laser driven 
flow and fluid manipulation, and for the stability of these. 



Appendix A: Equilibrium value after step-function switch-on 

We show here Eq. (asymptotic), the asymptotic equi- 
librium value at times r 3> 1 long after switch-on of a 
step-function laser, assuming no switch-off has occured. 



For r > 1, Eq. (33b) becomes 



P(£,T)-»2e-** / dxxe'* 2 \l (2&) - -Ji(2£ 



(Al) 



Using the following formulae from Ref. [37J 

^2 



Axx v+1 e' ax J„(/3x) 



P 1 



(2a) 



u+l 



-/3 2 /4a 



dxx 



v— 1 n — ax' 



J„({3x) =2»- 1 p-^(v,(3 2 /4a) 



where 7(0, x) = dte - *^ ^ 1 is the incomplete Gamma 
function, using J n (iz) = i n I n (z) and letting f3 = 2i£ we 
find 



J dxxe-* 2 1 (2£r) = -e« 2 , 



(A2) 



dxe~ x h{2£x) 



1 

2^ 



(1-e*). (A3) 



Inserted into ( Al ) this immediately gives Eq. ( 35 ) 
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